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The presence of a weak second-harmonic field in an intense-laser ionization experiment affects the 
momentum-resolved electron yield, depending on the relative phase between the uj and the 2uo component. The 
proposed two-color ’phase-of-the-phase spectroscopy’ quantifies for each final electron momentum a relative- 
phase contrast (RPC) and a ’phase of the phase’ (PP), describing how much and with which phase lag, respec¬ 
tively, the yield changes as function of relative phase. Experimental results for RPC and PP spectra for rare gas 
atoms and CO 2 are presented. The spectra demonstrate a rather universal structure that is analyzed with the help 
of a simple model based on electron trajectories, wave-packet spreading, and (multiple) rescattering. Details 
in the PP and RPC spectra are target sensitive and thus may be used to extract structural (or even dynamical) 
information with high accuracy. 

PACS numbers: 32.80.Rm, 34.80.Qb, 33.60.+q 


Momentum-resolved photoelectron spectra from strong- 
field ionization of atoms and molecules contain a wealth of 
information about the ionizing laser field, the target, and ul¬ 
trafast processes that may occur in it (see, e.g., Refs. mi 
for recent overviews). If one is able to disentangle this infor¬ 
mation one may use photoelectron spectra to image the entire 
ionization dynamics and the structural information convoluted 
into it. Thanks to the appealing possibility to analyze most 
of strong-field ionization dynamics in terms of semi-classical 
electron trajectories it has been demonstrated that many of the 
complex and, at first sight, puzzling features could finally be 
explained in simple and intuitive terms. Examples are various 
low-energy structures iMi, intra and inter-cycle interfer¬ 
ences Illll4l3h “holographic side-lobes” and the role of multi¬ 
ple returns |0-0], molecular strong-field ionization 117.18], 
or “interference carpets” [19]. 


Refined strong-field experiments make use of additional ex¬ 
perimental “knobs,” such as a pump-probe delay, a carrier- 
envelope phase (CEP) in the case of few-cycle pulses (see, 
e.g., Refs. 120] and 1 E 1 I E 2 I 1 for reviews), or a relative phase 
between two laser fields of different frequency |)23|-l29[]. The 
resulting changes in the photoelectron or, complementary, 
high-harmonics spectra may then help to unequivocally in- 
dentify certain ionization scenarios. In this Letter, we present 
experimental results obtained with a co-linear, two-color laser 
set-up. 


Consider the measurement of an observable (here the elec¬ 
tron yield). If this observable depends on a tunable parameter 
(here the relative phase between the c u and the 2cc pulse), a 
measurement of the change of the observable as function of 
the parameter reveals additional information about the under¬ 
lying physical mechanism (here the target and laser-sensitive 
ionization dynamics). A general question is how to repre¬ 
sent this additional information. If the parameter is periodic 
(like the relative phase) a Fourier-transform of the observable 
seems adequate. In our case, the momentum-resolved pho¬ 
toemission signal is Fourier-transformed with respect to the 
relative phase. The absolute value of the complex Fourier 
transform gives a relative-phase-contrast (RPC), its phase the 


’phase of the phase’ (PP). In this work, we present experi¬ 
mental RPC and PP spectra for various targets and analyze 
their common features in terms of “simple man’s theory” 
(SMT). The main findings are: (i) the overall structure of the 
PP spectra is largely target-independent and displays features 
that can be assigned to certain electron trajectories; (ii) target- 
dependent features are clearly visible in the RPC and PP spec¬ 
tra, thus making two-color PP spectroscopy an attractive ap¬ 
proach for revealing structural information. 

The basic structure of “ordinary” photoelectron spectra is 
qualitatively well understood. For linear laser polarization, 
the so-called “direct electrons” extend up to energies p^/2 = 
Aq/2 = 2 U p , where p z is the photoelectron momentum along 
the polarization direction, Ao is the vector potential ampli¬ 
tude, and Up = Aq /4 is the ponderomotive energy (atomic 
units are used unless otherwise stated). If electrons are driven 
back to their parent ion by the laser field, energies up to 10 U p 
may occur upon rescattering. A wealth of information is en¬ 
coded in this high-energy above-threshold ionization (HATI) 
part of the photoelectron spectra about both (i) the driving 
laser field and (ii) the structure of the target [3(1 |3JJ] (via its 
scattering cross section [32, 33]). Moreover, HATI is more ro¬ 
bustly accessible to theoretical modeling than the low-energy 
part of photoelectron spectra, which is plagued by the neces¬ 
sity to take Coulomb-corrections into account Il34l - [39ll . 


Bichromatic cc-2cc-pulses with parallel polarizations of the 
field components and adjustable relative phase <p were gen¬ 
erated with a set-up similar to the one in Ref. [ 40 ]. Briefly, 
a Tiisapphire laser system provides 100-fs pulses at 794 nm. 
The second harmonic is generated in a 100-/im thick BBO I 
crystal. The cc-2cc intensity ratio is controlled by detuning the 
phase matching conditions through a tilt of the crystal. Bire- 
fringent calcite crystals compensate for the time lag between 
the c 0 and 2cc pulse. The relative phase lag is controlled by 
two glass wedges mounted on piezo-driven motors. The laser 
pulses are focused into the extraction region of a homebuilt 
high-energy velocity-map-imaging (VMI) spectrometer 14111 
by a concave silver mirror having a focal length of 300 mm. 
Photoemission from Xe was used to optimize the pulse over- 
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FIG. 1: (color online). Schematic illustration of RPC P and PP <I>. 
Given a sequence of photoelectron momentum distributions (as ob¬ 
tained with the VMI spectrometer) for varied relative phase ip, the 
change in electron yield (dots) for given photoelectron momentum 
(small square in the spectra) is fitted by P cos(p + <f>) (red curve). 


lap. The MCP back plate of the detector system is gated using 
a fast electronic switch to suppress spurious signals. 

Our two-color field is described by the vector potential 

A(t) = Aoe z [smut + £sm(2ut + ip)] (1) 

in dipole approximation. The use of long pulses ensures that 
envelope effects are unimportant. Throughout this work the 
2ou field is kept weak (£ <C 1). 

FigureUintroduces schematically the two quantities—RPC 
P and PP —used to describe how the photoelectron yield 
Y (cp) changes as function of relative phase. For each final 
momentum p Z9 p x , the fundamental change in the yield, i.e., 
Yi(ip) = Pcos(ip + <f>), is extracted via the Fourier anal¬ 
ysis mentioned above. The fundamental (N = 1) compo¬ 
nent is found to dominate over higher-harmonic contributions 
Y/v> N > 1 for most final momenta. A similar analysis has 
been performed previously with respect to the CEP for SiC >2 
nanospheres ll42h . 

Experimentally determined RPC and PP spectra for Ar are 
presented in Fig. [2 The RPC in Fig. [2^ shows that the di¬ 
rect electrons vary most with ip (black and dark gray), while 
the rescattered electrons vary with about one order of magni¬ 
tude less contrast. The direct electrons with p z > 0 in the PP 
spectrum in Fig. [2t> behave predominantly sin ip to — cos ip- 
like (blue to black area labelled ’1’), the ones with Pz < 0 
behave cosip to — sin(^-like (green to red, ’2’). Most of the 
rescattered electrons behave phase-wise similar to the direct 
electrons in the opposite direction (’3’ red to black, and ’4’ 
blue to green). However, in the semi-circle-shaped momen¬ 
tum regions ’5’ and ’6’ beyond the respective 2 U p -cutoffs the 
electrons continue to follow phase-wise the direct electrons in 
the same direction, i.e., ’2’ and ’1’, respectively. This effect 
will be discussed in more detail below. Note that the overall 
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FIG. 2: (color online). Experimental RPC (a) and PP (b) spectra for 
Ar, calculated from 66 Abel-projected VMI spectra per p interval 
[0, 2tt\. The 794 nm / 397-nm components of the 100-fs two-color 
pulse had intensities ~ 10 14 /10 12 W/cm 2 , respectively (i.e., £ ®g 
0.05). In (a), P has been normalized to max P. Circles and vertical 
lines indicate 10 U p rescattering rings and 2U P cut-offs, respectively. 
White numbers in (b) are referred to in text. 


absolute phase is not determined by the experiment so that a 
cyclic shift of the color code for the experimental PP spectra 
is permissible. 


In Fig. [3] we show RPC and PP spectra for Kr, Xe, and 
randomly aligned CO 2 . Figures [ 2 J 3 and [34, f show that the 
overall structure of the PP spectra is universal: they resemble 
two overlapping clubs (indicated in Fig. [3), one colored in 
red to green (regions ’2’, ’5’, and ’4’ in Fig. [24, the other one 
blue to black (’1\ ’6’, and ’3’). The blunt parts of the clubs 
(regions ’3’ and ’4’) represent rescattered electrons. The tip 
regions ’5’ and ’6’ are investigated below. While the over¬ 
all structure is similar for all species, a target dependence is 
most obvious in the RPC spectra Figs.[3,c,e and[2k and also 
reflected in detailed features of the PP spectra. This allows 
two-color spectroscopy to be employed for imaging and as a 
sensitive test of theoretical models. In what follows, we aim 
at reproducing the results of Fig. [2] via simple modeling using 
SMT. 


SMT (see, e.g., 1221144[, [45D) should be able to reveal the 
origins of the common, overall features observed in the PP 
spectra. Given a vector potential (Q]) the electron momentum 
before rescattering reads p(t)/Aq =p(r) e z = a(r) — a{r\). 
The dimensionless time r = out , momentum p = p/Ao, and 
vector potential a = A/ Ao are introduced to highlight the 
universal scaling of SMT, and n = out \ with t\ the ionization 
time. Assuming isotropic scattering and considering the pla¬ 
nar motion in, e.g., the xz- plane, the momentum after reseat- 
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FIG. 3: (color online). Experimental RPCs (left) and PPs (right) for 
Kr (top), Xe (middle), and randomly aligned CO 2 (bottom). Laser 
intensities and color coding as in Fig. [2| The yellow dashed lines in 
(f) indicate the ’clubs’ discussed in the text. 


tering and once the pulse is off reads 

(£) = ('t 2 , ) + 'll)- ® 

Here, 72 = u;t 2 > t\ with 1 2 the rescattering time, and 6 
the scattering angle. The maximum momentum p max = y/5 
corresponds to the well-known cut-off energy 10 U p which 
occurs if ionization happens at at ti/ 27 t = 0.543 (or at 
0.043 in opposite direction) and the 6 = 7r-backscattering 
at T2/27T = 1.22 (0.72). The 10(7 P rescattering rings, i.e., 
rings in the momentum plane centered at —a(r 2 )e z with ra¬ 
dius |< 1 ( 72 ) — a(ri)|, are shown in Fig. [4^ below and indi¬ 
cated in all experimental and TDSE spectra. If electrons do 
not scatter during the first but the second return they end up 
in the opposite ^-direction and generate second-return rescat¬ 
tering rings. In particular, if ionization and rescattering occur 
at ti/ 27 t = 0.519 (0.019) and T 2 / 27 r = 1.74 (1.24), respec¬ 
tively, a local maximum cutoff energy of 7U P is reached, also 
indicated in Fig.|4^ (respective trajectory shown in Fig. 03). 

The fundamental club structure of the PP spectra is easy to 
reproduce within SMT: (i) looping over t\ and r 2 > t\ all 
electron trajectories that actually rescatter are considered, and 
(ii) weighted by W = W r W[ where W[ is an ionization and 
W r a rescattering probability. An instantaneous tunneling rate 
through a triangular barrier [46] W\ (ti ) ~ exp [—2/31 e (ti ) | ] 
(with e(r) = —d T a(r) the dimensionless electric field) has 
been taken. We are not interested in absolute numbers here, as 
our SMT modeling should be robust and depend only on a and 
r. In a first, crude approximation W v {j 2 — n) = (72 — r\)~ s 
with a spreading exponent s > 0 has been chosen. Such 
a rescattering probability accounts for wave-packet spread¬ 


ing (see, e.g., £L2,[44]) but neglects any momentum and an¬ 
gular dependence of the scattering cross section. If several 
trajectories end up in the same final momentum bin only the 
most probable has been considered, thus neglecting interfer¬ 
ence effects. The direct electrons have been weighted by 
Wf dir) (n) = Wi(n)exp(-/3^/|e(ri)|) with/? = 15 (to 
produce a reasonable lateral spread of the momentum distri¬ 
bution). The result of such a simple modeling for s = 3/2 
is shown in Fig.0:. The spreading exponent has been chosen 
smaller than s = 3 for free, 3D Gaussian wave-packet spread¬ 
ing in order to mimic Coulomb focusing [ 47]. The structure of 
the two overlapping clubs—one red, the other blue—is clearly 
visible. The weak (£ = 0.05) 2cc-component leads to a small 
extension of the rescattering cutoff (small strip of opposite 
color around the blunt club ends). 

Obvious disagreements between Fig.0: and all experimen¬ 
tal results are the oversimplified behavior consisting of es¬ 
sentially only two PPs (red or blue), and the too short low- 
momentum tip of the club ending at \p z \ = 1, corresponding 
to the 2(7 P cut-off for direct electrons. In all experimental and 
TDSE results the club tips extend beyond 2 U p (see regions ’5’ 
and ’6’ in Fig. [JJd). In order to reveal the origin of this mis¬ 
match we first note that in the SMT modeling of Fig. 0: the 
trajectories that rescatter at their first return dominate because 
of the (t 2 — Ti) -S -penalty for later returns. This is the reason 
why there is essentially only a sin-like or a — sin-like behavior 
for all final momenta (i.e., blue or red). In a refined modeling, 
a rescattering probability ~ W x (j 2 — Ti)cr(q, 0) with a differ¬ 
ential scattering cross section (in first Born approximation) of 
the form 


a(q,0) 


(2 Z) 2 

[fi 2 + 4 q 2 sin 2 ((9/2)] 2 ’ 


q = a(r 2 ) - a(ri) (3) 


for a screened potential V(r) = — Z exp(—fir)/r has been 
employed. The resulting spectrum in Fig. 01 is in much bet¬ 
ter agreement with the experiments. In fact, given the simple 
modeling, the agreement is striking. Going from Fig.0: to d, 
the PP in regions 3 and 4 changes from sin (blue) to — cos-like 
(black) for p z < 0, and from — sin (red) to cos-like (green) for 
p z > 0. The reason for this change lies in the introduction of 
an additional functional dependence on the relative phase via 
the cross section cr[q(ip), 0]. 

As in the SMT modeling all spectral features can be under¬ 
stood in terms of electron trajectories, we were able to identify 
the origin of the ring-like extensions ’5’ and ’6’ beyond 2(7 P . 
Figure [4^ suggests that these features are due to electrons that 
rescatter during the second return because the corresponding 
rescattering rings are centered around the 2U P cutoffs. A cross 
section of the form (0) favors forward scattering (i.e., small 
6), and the more so the larger the instantaneous scattering mo¬ 
mentum q = \cl{t 2 ) — a(r i)| is. As a consequence, the tra¬ 
jectories which scatter with lower-momentum later during the 
second return (incident from the opposite direction) may have 
a higher probability for large 9 than the trajectories that rescat¬ 
ter during the first return, despite the wave-packet-spreading 
penalty for late returns. It is thus the competition between the 
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FIG. 4: (color online), (a) Final SMT photoelectron momenta for 
£ = 0, absolute value of electric field at ionization time |e(n)| 
color-coded. Rings representing local maxima in final momentum 
| p z | due to rescattering during first, second and third return are la¬ 
belled. For better visibility, final momenta in respective opposite 
directions (originating from electrons emitted half a laser cycle ear¬ 
lier) are plotted light gray, (b) Normalized excursion a z of electron 
rescattering at first (solid) and second return (dotted), respectively, 
leading to final momenta indicated ’x’ in (a). SMT PP spectra for 
£ m 0.05, taking into account (c) ionization and wave packet spread¬ 
ing (s = 3/2) and (d) additionally a scattering cross section J3]) (with 
Z = IS and p =* 1/2). Same color coding in (c) and (d) as in previ¬ 
ous figures. White numbers in (d) analogously to Fig. 
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FIG. 5: (color online). RPC (a) and PP (b) spectra for Ar from TDSE 
(10 14 W/cm 2 , £ = 0.05). Durations of sin 2 -shaped pulses for 800 
and 400 nm were 60 and 50 fs, respectively. Additionally to 2 U p cut¬ 
offs and 10Up -rings, second-return cut-off rings are included, which 
form boundaries for features labelled ’A’ and ’B\ Slight left-right 
asymmetries visible in (a) are remnants of CEP-dependence. 


probability factors governing wave packet spreading and the 
momentum-dependent large-angle scattering which is respon¬ 
sible for the club tips beyond 2 U p . 

Figure \5\ shows RPC and PP spectra for Ar from a TDSE 
simulation 114311 in single-active electron approximation. Abel- 
projected photoelectron momentum spectra were calculated 
for different ip and subsequently treated like the experimen¬ 
tal spectra. The correct ionization potential of Ar was im¬ 
posed by choosing the 3po-orbital in the effective potential 
V(r ) = —[1 + 17exp(—2.11375r)]/r as the initial state. 
Comparison with Fig.[2]yields satisfactory agreement with re¬ 
spect to the overall club structure in the PP and the qualitative 
features in the RPC spectra. We note that the dynamic range of 
the experimental detection is two orders of magnitude lower 
than the one shown for the TDSE. This might be the reason 
why the marked features close to the second-return cut-off in 
the TDSE PP spectrum (labelled ’A’ and ’B’) are absent in 
the experiment. The remaining differences could be due to 
the idealized pulses and the neglect of focal averaging in the 
simulation, slightly different laser parameters in experiment 
and simulation, or the possibly inadequate (but computation¬ 
ally unavoidable) assumption of a single active electron in the 
TDSE calculations. In any case, such differences show the 
sensitivity of PP spectroscopy and its power to discriminate 
between various effects. 

In summary, we have introduced systematic two-color spec¬ 
troscopy based on RPC and PP spectra. The presence of the 
2cc-field component tags each emission time according to how 
ionization probability changes as function of relative phase. 


This change is subsequently mapped through the actual elec¬ 
tron dynamics to the final photoelectron momentum. Since 
several trajectories end up with the same final momentum the 
trajectory with the largest change in general dominates. We 
have revealed and explained the universal, overall structure of 
the PP spectra for various rare gas atoms and CO 2 . We further 
have shown that details in both the RPC and the PP are target- 
dependent, paving the way to employ them for imaging elec¬ 
tronic structure or dynamics, and as tests for theoretical pre¬ 
dictions or cross sections. We also anticipate that two-color 
spectroscopy will be able to discriminate among coherent, de¬ 
layed, or thermal electron emission, e.g., in non-sequential 
ionization of multi-electron systems, because the relative- 
phase dependence is maintained, affected, or destroyed, re¬ 
spectively. Any spectral feature identified in an “ordinary” 
photoelectron spectrum and hypothetically attributed to a cer¬ 
tain phenomenon (such as rescattering or focusing of certain 
trajectories, polarization of the target, or internal dynamics) 
can be further scrutinized, in that way testing whether the PP 
and RPC signatures of this feature are in accordance with the 
predictions from the conjectured theory. 

We thank M. Kling and S. Zherebtsov for bringing us in 
touch with the cc-2cc method and helpful support in the early 
stage of our studies. Fruitful discussions with M. Ivanov, com¬ 
puter time provided by the North-German Supercomputing 
Alliance (HLRN, project no. mvp00004), and support through 
the SFB 652 of the German Science Foundation (DFG) are 
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